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The field theoretic renormafization group is apphed to Kraichnan's model of a passive scalar 
quantity advected by the Gaussian velocity field with the pair correlation function oc 5{t — t') /k'^'^'^ . 
Inertial-range anomalous scaling for the structure functions and various pair correlators is established 
as a consequence of the existence in the corresponding operator product expansions of "dangerous" 
composite operators (powers of the local dissipation rate), whose negative critical dimensions de- 
■ termine anomalous exponents. The latter are calculated to order e"^ of the e expansion (three-loop 

' approximation). Submitted to Acta Physica Slovaca. 
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I In recent years, considerable progress has been achieved in the understanding of intermittency and anomalous 
scaling of fluid turbulence. The crucial role in these studies was played by a simple model of a passive scalar quantity 
advected by a random Gaussian field, white in time and self-similar in space, the so-called Kraichnan's rapid-change 
There, for the first time the existence of anomalous scaling was established on the basis of a microscopic 
and the corresponding anomalous exponents were calculated within controlled approximations [3-6] and a 
systematic perturbation expansion in a formal small parameter |0] . Detailed review of the recent theoretical research 
on the passive scalar problem and the bibliography can be found in ||^. 
^ The advection of a passive scalar field 6{x) = 9{t,x.) is described by the stochastic equation 

d ■ dt0 + M)9 = i^oAe + /, (1) 

(~| , where dt = d/dt, di = d/dxi, vq is the molecular diffusivity coefficient, A is the Laplace operator, vix) = {vi{x)} is 
the transverse (owing to the incompressibility) velocity field, and / = f{x) is an artificial Gaussian scalar noise with 
, , zero mean and correlator 

{f{x)fix'))=Sit-t')C{r/L), r=|x-x'|. (2) 

in ■ 

The parameter L is an integral scale related to the scalar noise and C{r/L) is some function finite as L oo. Without 
I/"} \ loss of generality, one can take L = oo and C(0) = 1. 

• In the real problem, the field v(a;) satisfies the Navier-Stokes equation. In the rapid-change model it obeys a 
Gaussian distribution with zero mean and correlator 

O : 

{v^{x)vj{x')) = Z?o ^1^^ j rfkPy(k) k^''-' exp[ik • (x - x')], (3) 
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where -Pij(k) — Sij—kikj/k^ is the transverse projector, fc = |k|, > is an amplitude factor, d is the dimensionality 
of the X space and < £ < 2 is a parameter with the real ( "Kolmogorov" ) value e = 4/3. The infrared (IR) 
regularization is provided by the cut-off in the integral from below at fc ~ m, where m = 1/f is the reciprocal 
$H of another integral scale £; the precise form of the cut-off is unessential. The relation Dq/vq = A^ defines the 
characteristic ultraviolet (UV) momentum scale A. 

The issue of interest is, in particular, the behaviour of the equal-time structure functions S'„(r) — ([0(i,x) — 
6'(i,x')]"), r = |x — x'l, in the inertial range A ^ 1/r ^ m. (With our assumptions, the odd structure functions 
vanish; they become nontrivial if the correlation function (vf) is nonzero or a constant gradient of the scalar field 
is imposed.) If we neglect the advection, that is, the nonlinearity in (|l|), these functions can be easily calculated: 
Snif) — Cn^Q^r^" with some constants C„. In the full nonlinear problem, they become dependent on two additional 
variables Ar and mr: 

5„(r) = i^o""^'"^«(Ar,mr). (4) 

The analysis by Refs. shows that the behaviour of the functions F„ in (^ at Ar ^ oo, mr (inertial range) 
is nontrivial: 

Sn{r) ~ C„i/o-"r2"(Ar)-""(mr)~'3", (5) 
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that is, the dependence on the both scales persists. Moreover, the exponent /?„ is a nonhnear function of n, the 
phenomenon referred to as "anomalous scaling" in the theory of turbulence. 

Relations (|^) and explicit expressions for the exponents and /3„ up to the first order in 1 / d and e were derived in 
[^,||. Within the "zero- mode approach" to Kraichnan's model, developed in those papers, anomalous exponents are 
related to the zero modes (unforced solutions) of the closed exact equations satisfied by the equal-time correlations. 

In Ref. 01 and subsequent papers [9-17], the field theoretic renormalization group (RG) and operator product 
expansion (OPE) were applied to the rapid-change model and its descendants. In the RG approach, anomalous 
scaling emerges as a consequence of the existence in the model of composite operators with negative scaling dimensions, 
identified with the anomalous exponents. This allows one to construct a systematic perturbation expansion for the 
anomalous exponents, analogous to the famous e expansion in the RG theory of critical behaviour, and to calculate 
the exponents to the second ||^,^,^ and third orders. For passively advected vector fields, any calculation 

of the exponents for higher-order correlations calls for the RG techniques already in the lowest-order approximation 
[15-17]. Besides the calculational efficiency, an important advantage of the RG approach is its relative universality: 
it is not related to the aforementioned "solvability" of the rapid-change model and can also be applied to the case of 
finite correlation time or non-Gaussian advecting field [^ . 

Below we briefiy review the main ideas and results of the RG approach to the rapid-change model; detailed exposition 
and more references can be found in l^, [l^ , |l^ , |l7| . 

The solution proceeds in two stages. In the first stage, the RG equation is employed to find the asymptotic behaviour 
of the functions in (^) with respect to their first argument (Ar — > cxi) at fixed mr. This gives 

Sn{r) ~ z.o-"r2"(Ar)-""F„(mr). (6) 

The total power of r, namely 2n — an, is (up to the minus sign) the critical dimension of the quantity on the left-hand 
side. However, the form of the function F„(mr) cannot be determined by the RG equation. 

Expression (^ is the special example of the general property established by the RG method for problem (0)-(|^): 
the property of the IR scale invariance. Similar asymptotic expressions can be obtained for all correlation functions; in 
particular, for the equal-time pair correlation functions S'„fc(r) = ("I>„(a;)"I>fe(a;')) of the monomials $„ ~ [di9{x)di9{x)\'^ 
one obtains: 

5„fc(r) ~ v-^-\Kr)-^--F^k{mr). (7) 

In the RG technique, the exponents a„ and ^nk in (^, (0) arc calculated in the form of series in the parameter e from 
(||), which is therefore the analogue of the quantity e = 4 — d in the (f/^ model of critical behaviour. 

From the general RG viewpoint, model (|l|)-(^ is simpler than the conventional 0"* model in a few respects: the 
coordinate of the fixed point is exactly determined by the one-loop approximation, renormalization of the primary 
field 9{x) and its powers 0"(x) ("composite operators" in the field theoretic language) is absent; the "mass" m is not 
renormalized. Therefore the exponents a„ in (j^) are found exactly, a„ — ne, that is, they have no corrections of order 
e^, £^ and so on. However the exponents 7„fc, given by the sums 7„fc — A„ -I- of the critical dimensions A„ of 
the composite operators $„, have nontrivial e expansion. The knowledge of these dimensions allows one to find the 
asymptotic form of the mean values of the $„ at m/A 0: 

($„) -a.;^o^"(m/A)^". (8) 

Relation (H) gives a contribution in ($„) that is nonanalytic in m/A; this contribution prevails if A„ < (the 
operator is "dangerous"), which is indeed the case: in the first order in e one obtains A„ = —2n{n — \)£/{d + 2), 
cf. [^. Strictly speaking, because of the mixing in the renormalization the monomial becomes a finite linear 
combination of the "scaling" operators with definite dimensions A^ with k < n; in representations (Q) and (^ we 
show only the leading term with the maximum |Afe|. 

The scaling functions Fn{mr), i^„fc(mr) in (|^), (^ can be calculated as series in e, but such representations become 
useless if the behaviour of these functions at mr — > is studied: the actual expansion parameter then becomes 
£ln(mr) rather than e. Summation of such "large IR logarithms," in contrast with the "large UV logarithms" with 
In(Ar), lies beyond the scope of the RG method in the narrow sense of the word, because the forms of the scaling 
functions are not determined by the plain RG equations. Like in the theory of critical behaviour, the behaviour at 
mr —> can be extracted from the general solution of the RG equations by means of the operator product expansion 
(OPE). According to the OPE, the product of two (renormahzed) composite operators Oi(xi)02(x2) (for example, 
6" or renormalized analogues of the monomials $„) a.t ti = t2 = t, r = xi — X2 ^ and the fixed "centre of mass" 
X = (xi + X2)/2 — const is represented in the form 

Oi(a;i)02(x2) =^C„(r)0„(x,t), (9) 
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where Cq are coefRcients analytic in (mr)'^ and the summation, in general case, runs over all local composite operators 
Oa{x) with definite critical dimensions A^. Correlation functions like (|^), (0) are obtained by averaging relations 

the mean values {Oa{x)) oc m^" appear on the right-hand sides, they give rise to the contributions of the form 
(mr)'^'' in the scaling functions Fn{mr), EnkiTtir). The leading terms at mr — > are determined by the operator 
with the minimum possible value of Aq,. In conventional models like 0^ the leading term is related to the simplest 
operator Oa{x) = 1 with A^ = 0. However, in model the critical dimensions of all operators $„ are negative 

(see above), and they determine the leading terms for mr ~f 0. Although the operators containing fields without 
the derivative d also have negative dimensions, they do not contribute to the OPE for the quantities like (^, ^ due 
to the invariance of the latter with respect to the shift 0^6 + const of the field 0{x). 

If the operator product expansion for any given function or (0) included contributions from all the operators 

we would have to sum all of them in order to obtain the behaviour for mr — > 0, because the spectrum of their 
dimensions is not bounded from below (the most dangerous operator does not exist). This is the case in more realistic 
models of fully developed turbulence; see [18-20]. Fortunately, this problem does not exist in model (|l])-(||) due to 
the linearity of equation (|^) in 0{x): one can show that the number of fields in any composite operator Oa{x) that 
enters the right-hand side of representation (||) cannot exceed the total number of fields on the left-hand side. We 
thus arrive at the following expressions for quantities (^, (j^): 

Sn{r) ^ i.o-"r2"(Ar)-"^ ^ C„,,(mr)^% (10) 



5„fe(r) ~ i.o-"-'=(Ar)-^"-^'= J2 CnkAmr)'-' (H) 

with coefficients C„^s, Cnk,s dependent on e and d and corrections of the form (mr)^+'^^'^\ Since the dimensions A^ 
are negative and decrease as s increases, the leading asymptotic term in (|l^, ( pT| ) is given by the contribution with 
the maximum s, that is, with s = n in (|lo| ) and s — n -\- k in (|Tl|): 

5„(r) ~ C„z.o-"r2"(Ar)-"^(mr)^", (12) 



Snk{r) ~ C„fcZ.o""-''(Ar)-^"-^H"^0^"^'=- (13) 

The critical dimensions are calculated as series of the form A„ = X^fc^i £'^Ai''''jthe first-order term was already 
given above. We have calculated the dimensions A„ in the second 0] and third [|3|jlJ] orders of the e expansion. The 
second-order result has the form: 



A(2) = 



n{n — 1) 



-A{d+\){d + Af + 



(d-l)(d + 2)3(d + 4)2 

-t-3(d - l)(«i + 2)(d + 4)(d + 2n)h{d) - 4(rf + l){d + 2){d+'in- 2)h{d + 2)1 , (14) 

where h{d) = F{1, 1; d/2 + 2; 1/4) and F{- ■ ■) is the hypergeometric series. Simpler expressions are obtained for integer 
d, in particular, h{2) = 8[1 — 31n(4/3)], h{3) = 10{n\/3 — 16/3), while for the other integer d analogous expressions 
can be obtained from the recurrent relation 3h(d) + {d+2)h{d+2)/{d+4:) — 4. The third-order coefhcient is presented 
in |l|jl|. 

— — (2) 

From the above formulas one obtains Aj^ ^ = in agreement with the exact result Ai — 0. The formal proof of the 
latter is based on certain Schwinger equation, which has the meaning of the conservation law for the "energy" 9^{x); 
see 0. This means that the second-order structure function S2 is not anomalous in agreement with the well-known 
exact solution obtained in 

We note that the family of operators is "closed with respect to the fusion" in the sense that the leading term of 
the correlation function ($n^m) is given by the operator ^n+m from the same family with the summed index n + m. 
This fact along with the inequality A„ -t- A,„ > An+rm which is obvious from the explicit expressions for A„, can be 
interpreted as the statement that the correlations of the local dissipation rate in the model (|l|)-(^ exhibit multifractal 
behaviour, cf. ^j^. 

An important issue that can be discussed on the example of the rapid-change model is that of the nature and 
convergence of s series in models of turbulence and the possibility of their extrapolation to finite values of e ^ 1. 
The knowledge of the three terms of the e expansion in model (|^)-(^) allows one to discuss its convergence properties 
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and to obtain improved predictions for finite e in reasonable agreement with the existing nonperturbative results: 
analytical and numerical solutions of the zero-mode equations and numerical experiments | 

The RG and OPE approach presented above can be generalized to the cases where compressibility, anisotropy or 
finite correlation time are present [10-14] and the passive advection of vector (e.g., magnetic) fields [15-17]. 

Let us conclude with a brief discussion of the RG and OPE approach to a more realistic model of fiuid turbulence: 
the stirred Navier-Stokes equation; see [18-20]. Dangerous operators in that model are absent in the e expansions 
and can appear only at finite values of e. This means that they can be reliably identified only if their dimensions 
are derived exactly with the aid of Schwinger equations or Galilean symmetry. Due to the nonlinear nature of the 
problem, they enter the corresponding OPE's as infinite families whose spectra of dimensions are not bounded from 
below, and in order to find the small-mr behavior one has to sum up all their contributions in the representations 
like (p^, (p^). The needed summation of the most singular contributions, related to the operators of the form w" 
with known dimensions, was performed in [ p^ using the so-called infrared perturbation theory for the case of the 
different-time pair correlation functions. It has revealed their strong dependence on m, which physically can be 
explained by the well-known "sweeping effects." This demonstrates that, contrary to the existing opinion [ p5| , p6t , 
the sweeping effects can be properly described within the RG approach, but one should combine the RG and OPE 
techniques and go beyond the plain e expansions. Analysis of the m dependence of the Galilean invaiant objects like 
the structure functions requires the explicit construction of all dangerous invariant scalar operators, exact calculation 
of their critical dimensions, and summation of their contributions in the corresponding OPE. This is clearly not a 
simple problem and it requires considerable improvement of the present techniques. 
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